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Abstract. In this note we discuss the relationship between the 
generating functions of some Hermite polynomials H , 
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generalized Airy-Heat equations 
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^- J exp ja(iA)" - ^X 2 t + iXx^ dX, 

higher order PDE's 
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and generating functions of higher order Hermite polynomials H ^ : 
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In particular, we show that under some conditions, these problems 
are equivalent. 



1. Introduction 

This note is motivated by the observation, made in Hernandez- del- 
Valle (2010), that all the known closed form densities of hitting times 
of Brownain motion may be expressed in terms of heat polynomials, 
which may be alternatively formulated as Hermite polynomials. 

Except in the case of the constant and linear boundaries, it seems 
that solutions of the hitting time problem will involve infinite series of 
Hermite polynomials. In general, these series expansions will not be 
convergent, as has been pointed out in Widder and Rosenbloom (1959). 

Hence, in the present work we derive an alternative representation, 
which relates to generalized Airy functions, higher order Hermite poly- 
nomials, and higher order PDE's. 

l 
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2. Triple Lacunary generating function of Hermite 
polynomials and alry-heat equation 

In Gessel and Jayawant (2005) the authors derive — using umbral and 
combinatorial arguments — a generating function for Hermite polyno- 
mials of order 3, namely: 

(i) IJ 1 ^)^ = -jttmT 2Fo U' e ; " : (i - swzy ) 



in which w — (1 — a/1 — 12wz)/6z = u ■ C(3uz), where C(x) — (1 — 
x is the Catalan number generating function, and 2-P0 is the 
so-called 'mysterious' hypergeometric function, which in turn is related 
to the Bessel function B of order 1/3 through 



^o(^;-;e) = ^e-^ 1/3 \~ 



and alternatively, By^ is (under certain conditions) equivalent to the 
so-called Airy function Ai, defined as: 

(2) Ai(x) := — / exp {iX + iXx} dX, 

2tt J_ qo 

where i := \/— 1. 

2.1. Alternative derivation of (0Q). A solution to the forward heat 
equation 

<9iu , . 1 d 2 w , . 
^' x) = 2^ ( *' x) 

is given in terms of the following inverse Fourier transform: 

1 r+00 3 

(3) w ( t ,x) = — e^-^+^dX. 

2vr J_ 00 

After completing the cube in the exponential function, equation ([3]) 
may be written as: 

ft 3 tx 1 / t 2 
iu(t, x)=exp|- + -|Ai (a; + - 

where Ai is defined in (TSl). 
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Finally, equations (EE]) and ([3]) are 'symbolically' equivalent, by inter- 
changing the McLaurin expansion of e lA3//3 in ([3]) 

i r+oo 3 

(4) w(t,x) = — / e i2 T^ x2t+iXx d\ 
2tt 
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(5) = £ 



-3) 



-j 
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and observing that the term within the squared brackets correspond 
to the so-called derived heat polynomials of order n — introduced by 
Widder and Rosenbloom (1959): 



1 



+oo 



cu n (t,x) := ^- / e lXx -^ x2t (-i\) n d\ 



OO 



(6) = t-^k(t,x)H n ^y 

It follows that (jSJ), together with (j7]) is equivalent to 



(7) W(t,x) = V^3n(«)^. 

j=0 



As a corollary, we conclude that generating functions of Hermite 
polynomials of the form 

OO j 

are equivalent to Airy-heat functions, defined as: 

— I exp(aA m --\H + i\x\ dX, 

with the appropriate identification of variables u and z. 

There is though an important issue, the interchange of summation 
and integration in ([5]) is only symbolic, since e tXy / 3 — although entire — 
is not of order (2, a). In other words, approximating (jlj) with ([5]) is not 
allowed. 
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3. Correspondence between expansions of Hermite and 
Higher order Hermite polynomials 

In this section we describe a sort of 'duality' between expansions of 
the form: 



3=0 

and generating functions of higher order Hermite polynomials H^: 

OO j 

3=0 J ' 

defined as: 



Hf\t,x) = e td "x j 



= ^77 



, U ~nk)\kV 

But first, let us embed solutions of the heat equation into higher 
order PDE's: 

3.1. Higher order PDEs and the heat equation. Let us start by 
defining the following bivariate PDE 

. . du d n u . . 1 d 2 u . . 

(8) m^ x) = a ^ x) + r^ x) 

where a and s are arbitrary coefficients. It follows immediately that it 
has particular solutions such as: 

u(t, x) = exp | ^a\ n + ^s\ 2 ^j t + \x 

= exp | (^a(i\) n - 7^ s ^j t + i\x 

= exp | ^a(i\) n - ^sX 2 ^ * j sin(A:r) 

= exp | ^a(iX) n - ^ s ^j *| cos(Ax) 

or 

(9) u(t,x) = e xt v(x) 
where v is a solution to the following ODE 

\ v (x) = av {n \x) + \sv {2 \x). 



AIRY-HEAT, HERMITE AND HIGHER ORDER HERMITE 5 

Which follows after substitution of ([9]) in (JSJ). 

Alternatively, applying the Fourier transform with respect to the 
space variable x to ([8]): 

u t = a{iX) n u — 2 S ^ 2 ^ 

(ai\) n - ^sA 2 a. 

and solving for t, leads to 

1 



u(t, A) = u(0, A) exp ia(iX) n t - -sX 2 t 

Hence, in terms of the inverse Fourier transform equation (JSJ) admits 
solutions of the following form 

1 f +co ( 1 1 

(10) u(t, x) = — J w(0, A) exp |a(iA) n t - -sA 2 t + zAxj dA. 

For example, in solution ( fTOl) let a = — 1/3 and u(0, A) = 1 

u(t, x) = — / exp < i sX t + iXx > 

which after completing the cube is equivalent to 

(11) „(t, x) = exp + f } ^ 1/3 Ai (V 1 / 3 + x}) • 

In particular, when t — 1, then ffTTl) is a solution to the heat equation 
in variables and alternatively when s = 0, it leads to an Airy 

kernel, introduced by Widder (1975). 

Note that, either by a proper time change or simply by setting t = 1 
in CEO} 

u(l,x) = — [ u(0, A) exp \a(iX) n - -sA 2 + iXxX dX 

leads to solutions of the heat equation in terms of variables (s,x): 

dv . . 1 d 2 v . . 

On the other hand if we set s = and u(0, A) = 1 in ( jTUl) 

^ r+oc 

i(t, x) — — / exp {a(iX) n t + iXx} dX 
2tt 
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we get solutions to higher order PDE's which may be alternatively 
described in terms of higher order 'Airy' functions: 

AiW (u) := — / exp {a(iA) 3+n + tXx] dX 
2tt > /_ 00 

(that is, when n = we obtain the classic Airy function), i.e. 



( 12 ) = tf&S)*™ 



x 



t l/(n+3) J • 



In general, with u(0, A) = 1, it follows that fTiTJj) is a convolution of 
T2l) and the heat kernel: 



/ N 1 f ^ 

(13) vis, x) = . exp < 

We may now describe the duality between generating functions of Her- 
mite and higher order Hermite polynomials. 

3.2. Duality between generating functions of Hermite and higher 
order Hermite polynomials. As described in the previous subsec- 
tion a solution to 

du . . d n u . . 1 <9 2 u . . 

is given by 

1 /" +00 f 1 1 

x) = — / exp < a{i\) n t s\ 2 t + iXx > dX. 

2tt y_ 00 I 2 J 

which can be alternatively seen as the convolution 

p+oo 

(14) u(t,x):= v(T,y)w(t,x - y)dy 



oo 



between the fundamental solutions of the following pair of PDE's: 

dw d n w 

[see equations ( IT2l) and ( fl~3l) respectively] where we have set a = 1 and 

T = t ■ S. 

On the other hand recall the definition of the higher-order Hermite 
polynomials: 

(16) H^^x) = e t9 "x j 
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which after differentiation with respect to t on both sides of (ITBj) we 
have that 

d t H<f\t,x) = d n x [e t9 "x j ] 



is a solution of f[T5]) and at time t = 

Hj(n)(0, x) = x 

i.e. it solves a Cauchy problem. After substiution of v equation f[T51) 
in (fH| and subsequently taking its MacLaurin expansion, it follows: 

[ +co 1 f Z/ 2 1 

u(t,x) := / -j=eKpl-—\w(t,x-y)dy 



-oo 

1 00 / I \ ./ I + N 



J! 



i ^ /_j_y i ( } 



y 2j w(t,x- y)dy 



j=0 

So for instance, a symbolically equivalent triple lacunary generating 
function described in ([7]) based upon higher order Hermite polynomials 
is given by: 

^f:(--Y^^x). 
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